General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



NASA TECHNICAL nasa tm 78,434 

MEMORANDUM 

(NASA-TM-78434) EfPECIS Of UNSTEADY N77-j2078 

ASBOEYNAWICS CN fCTOo AEBOEiASTlC STABILITY 
(NASA) 30 p fcC A02/MP A01 CSCL 01A 

UiiClas 

G3/02 47866 


» 

t'. 

S 


< 

CO 

< 

Z 


EFFECTS OF UNSTEADY AERODYNAMICS ON ROTOR 
AEROELASTIC STABILITY 


D(>nald L. Run/ 


Ames Research Center. NASA 
and 

Ames Directorate. USAAMRDL. AVRADCOM 
Ames Research Center. Moffett Field, Calif. 94035 


4 


September 1977 


-V C S 5 
- .. 

'c ■ 



'■iLirf 




■'c,i-ri7V 



EFFECTS OF UNSTEADY AERODYNAMICS 
ON ROTOR AEROELASTIC STABILITY 

Donald L. Kunz 

Research Scientist 

Ames Research Center 
Ames Directorate 

U. S, Ainny Air Mobility R&D Laboratory 
AVRADCOM 

Moffett Field, California 
Abstract 

An analysis was conducted to study the effects of unsteady aerodynamics 
on the stability characteristics of helicopter rotor blades. A simple physical 
model of each blade was used together with Theodorsen, Loewy, and quasi-steady 
aerodynamics to derive the equations of motion. The stability analysis com- 
paring the effects of using each of the three theories revealed some signifi- 
cant differences between the Loewy and Theodorsen results. These included 
increases and decreases in lead-lag damping, localized around integer lead-lag 
frequencies. It was also shown that the standard method of multi-blade coordi- 
nates must be modified for use in conjunction with Loewy aerodynamics. 
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Linear, two-dimensional lift curve slope 
Tip loss factor 

Theodorsen's lift deficiency function 
Loewy's lift deficiency function 
Profile drag coefficient 
Blade chord 

Nondimens ional blade chord, c/R 
Airfoil section drag force 

Hinge offset nondimens ionalized by R, Pig. 1 
Wake spacing nondimensionalized by R, Pig. 2 
Blade inertia about the hinge 

Blade torsional inertia nondimensionalized by I 

Blade geometry correction factors, Eq. (16) 

Combined flap and lead-lag spring stiffness at 
9^ = 0. = 0, Eqs. (AlO-All) 

b Cl 

Plap and lead- lag spring stiffnesses at the blade root. 

Pig. 1 

Plap and lead- lag spring stiffness of the inclinable hub 

springs , Pig . 1 

Reduced frequency, (iJc/2E2r 

Airfoil section circulatory lift force 

Airfoil section noncirculatory lift force 

Nondimens ional aerodynamic hinge moments 

Nondimensional elastic hinge moments 

Nondimens ional inertial hinge moments 
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Modal frequency ratio , o)/fl 
Number of blades 

Uncoupled rotating flap frequency 
Generalized lift deficiency function 
Integrated lift deficiency functions, Eq. (17) 

Blade number (0, 1, 2, . . . ) 

Rotor radius 

Elastic coupling parameters, Eqs. (A8-A9) 

Distance outboard of the hub centerline nondimensionalized 
by R 

Modal eigenvalue, b + iw 

Blade velocity in the plane of rotation, Eq. (12) 

Distance outboard of the hinge nondimensionalized by R 
Blade root cutout (from the hinge) nondimensionalized by R 
Blade velocity normal to the plane of rotation, Eq. (11) 
Blade flapping deflection. Fig. 1 
Lock number, pack'^/l 
Stiffness parameter, Eq. (7) 

Blade pitch angle with respect to V 
Blade lead- lag deflection. Fig. 1 
Blade geometric pitch angle. Fig. 1 
Inclinations of the blade and hub principal flexural 
axes , Fig . 1 

Pitch-flap coupling parameters, Eqs. (20-22) 

Pitch-lag coupling parameters, Eqs. (20-22) 

Air density 
Modal damping 
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Lead-lag damping 
Inflow angle 
Interblade phase angle 
Rotor rotational speed 
Modal frequency 
Lead-lag frequency 

Nondimens ional, uncoupled flap, lead-lag, and torsion 
natural frequencies, Eqs. (A12-A13) 

Equilibrium and perturbation quantities 
[d( )/dt]/fi 
[d^( )/dt^]/n^ 
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Introduction 


Over the past few years, many investigators have studied the aeroelastic 
stability of hingeless helicopter rotor blades. For the most part, these 
researchers have used quasi-steady, strip theory aerodynamics to develop the 
aerodynamic forces [1-3]. Although the results obtained from such analyses 
have often correlated well with experiment [A]j and have provided valuable 
Insights into the physical problem, some discrepancies have appeared which 
cannot be explained within the limits of the theory [5]. A possible source 
of these discrepancies is unsteady aerodynamics. 

The aerodynamic theories used in this investigation include Theodorsen's 
unsteady theory [6], Loewy's rotary-wing theory [7], and a quasi-steady approxi- 
mation to the two unsteady theories. By comparing the results obtained by 
using each of these aerodynamic formulations, the relative effects of unsteady 
flow and the shed wakes beneath the rotor plane can be determined. 

While other investigators have used unsteady aerodynamics in stability 
analyses [8-10], their main concern has generally been the determination of 
stability boundaries. The primary reason for the limited scope of these 
analyses has been that the aerodynamic theories only apply for simple harmonic 
motion. It is the intention of this study to treat transient motion, by assum- 
ing that motion to be nearly simple harmonic, as well as the case of neutral 
stability. In addition, the work presented here observes the effects of the 
shed wakes below the rotor plane on the behavior of a multibladed rotor. 

Blade Model 

The physical model chosen to represent each individual blade of the rotor 
was purposely kept simple, so that the aerodynamic effects would be obscured as 
little as possible by the blade dynamics. It is the same as the model used for 
the theoretical development in Reference 5 . 
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As shown in Figure 1, the blade model has flap and lead-lag degrees of 
freedom, perpendicular and parallel to the plane of rotation. While the blade 
is rigid in torsion, provisions are made for including kinematic pitch-lag, 
pitch-flap, and elastic coupling. Hinge offset can also be represented. In 
addition, the hinge sequence can be changed from lag-flap to flap-lag by 
rotating 9^ from 0“ to 90“ . All of the blades are attached to a fixed hub 
which rotates at a constant speed, fl. 

Aerodynamic Tb r nries 

Before proceeding with the derivation of the equations of motion, the 
aerodynamic theories that are used will be briefly discussed. Since both the 
Theodorsen and Loewy theories are widely known, this discussi?-n will serve only 
to compare their mathematical models and approaches . It will also point out 
the limitations their use imposes on the stability analysis. 

The unsteady aerodynamic theory developed by Theodorsen [6] considers a 
thin airfoil undergoing simple harmonic motion, and being trailed by a straight 
wake. The vortex strength of this wake varies sinusoidally with times, and 
results in a lift deficiency function, C(k) , which is associated only with the 
circulatory portion of the lift. This lift deficiency function is only depen- 
dent on the reduced frequency of the airfoil oscillations . 

Loewy [7] modified the wake model in Theodorsen 's analysis to more closely 
represent the casu of a rotating wing. As shown in Figure 2, the airfoil 
section at any spanwise station is trailed by a wake in the plane of the rotor, 
as in Theodorsen’ s model. However, Loewy also includes vortex sheets below the 
plane of rotor, which result from the vorticity shed by previous passes of the 
blades being carried downward by the induced velocity. For ease of computation, 
the wake sheets are assumed to be parallel to the plane of rotation and without 
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Co the plane of rotation and without curvature. They are also assumed to 
extend infinitely far ahead of and behind the airfoil, as well as below it. 


Loewy's analysis of this aerodynamic model follows the same approach as 
that of Theodorsen, and the results are similar. The lift deficiency function 
C(k,ra,h,N,i|^^) , is again only associated with the circulatory lift. However, it 
is now dependent on the frequency ratio, the wake spacing, the number of blades, 
and the interblade phasing, as well as on the reduced frequency. 

The quasi-steady approximation to the Theodorsen and Loewy theories is 
obtained by letting the lift deficiency function equal unity. This has no 
effect on the non-circulatory lift, but does make the circulatory lift indepen- 


dent of blade frequencies and wake geometry. 

Since the lift deficiency functions of the unsteady theories are dependent 
on reduced frequency, a constraint is placed on the applicability of the sta- 
bility analysis. That is, only those modes which have frequencies corresponding 
to the reduced frequency used to compute the value of the lift deficiency func- 
tion are valid results. Therefore, unlike when the quasi-steady theory is used, 
a separate calculation must be made for each mode. 

Another limitation on the validity of the stability analysis is applied 
when one is considering transient motion. Both Theodorsen's and Loewy's 
theories were developed only for simple harmonic airfoil motion. However, for 
lightly damped modes, the airfoil motion will be nearly simple harmonic. 


gSt ^ ^(a+i(a)t ^ (1) 

e*^ = 1 + cf + a^/2! + . . . + 0 '^^/n! 

Cf s t 

Therefore, when the damping is small, e is approximately unity and e is 
approximately e^^^. This assumption of nearly simple harmonic motion 
essentially means that only the results for the lead-lag mode are valid, since 
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the damping of the flap mode is generally outside the limits of this 
approximation . 

One final problem encountered when dealing with the Theodorsen and 
Loewy lift deficiency functions results from their being represented by 
complex numbers. If the resultant complex aerodynamic coefficients are not 
expressed in real form, the two eigenvalues for each mode will not necessarily 
be complex conjugates. The problems in physically interpreting such a solution 
are obvious. It is therefore necessary to invoke the assi'mption of simple 
harmonic or, in this case, nearly simple harmonic motion in order to express 
the aerodynamic coefficients as real numbers. 


Equations of Motion 

From Reference 5, the inertial and elastic moments about the hinges are 


where 
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Since all of the aerodynamic theories used in this study are similar in 
form, the aerodynamic moments are derived for an arbitrary lift deficiency 
function, Q. The expressions for lift and drag are based on Reference 11, 
with the small pitch angle assumption removed. 


nc 


pac 
2 4 


z cos 0 + ze sin 9 + V sin 6 + Ve cos 6 + "^ £ 




L = V I z cos 0 + V sin e + 
c z 


f^} 


B . £|S. v2|_a 

2 \ a / 


where 


z = - j^i'J(|)^(e + x) + 3 x^R 
V = I^J^Ce + x) + ^x^R 


e = sin 3 


( 8 ) 

(9) 

( 10 ) 


( 11 ) 

( 12 ) 

(13) 


Note that the blade vertical velocity, z, contains contributions from the 
induced velocity and blade flapping; while V, the horizontal velocity, has 
contributions from the rotor rotation and the blade lead- lag motion. The 
terms containing £ are a result of the virtual pitch rate induced by the 
coning of the blade as it rotates. 

To obtain the aerodynamic moments. Equations 11-13 are substituted into 
Equations 8-10, and the lift and drag resolved along the flap and lead-lag 
axes. The forces at each blade section are then multiplied by their moment 
arms about the hinge, and the product integrated over the span of the blade. 
The resulting moments are 

fig = ^ I (^2®^ ^ K^) (sin 0 - cos 0) - cos oj B 

+ || sin 2e'| Z - [(Q2e - Q3) cos 0] B 
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+ {(^2^ Q 3)(2 sin 0 - (f)^ cos 0 )^ 


2(Q,e + Q^) 


+ (K^e + K2 )cos 0 


( 14 ) 


M (K2e +2K3e + K^) 
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(15) 


where 


K = — (B-.e)'^ - x^l 

n n L oj 


(16) 


B-e 

Q = 4 / x^Qdx 

n J 

X 

o 


(17) 


The flap equation is obtained by summing Equations 3, 5, and 14, and 
setting the result equal to zero. Equations 4, 6 , and 15 are summed and set 
equal to zero in order to obtain the lead-lag equation. By perturbing the 
blade about an equilibrium position, the equilibrium and perturbation 
equations are obtained. 


3 = 3 + A3 

0 

(18) 

C = Cq + A? 

(19) 

0 = 0 + 0qA3 + 0,-AS 

op ^ 

( 20 ) 

'b"®b 

( 21 ) 
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( 22 ) 


0, = 0, + 0_ Ag + 0„ AC 


The complete equations can be found in the appendix. 


Multl-Bladed Rotor 

When working with rotors having more than one blade, the standard practice 
is to transform the isolated blade equations into multi-blade coordinates [12]. 
Under this transformation, a four-bladed rotor would have four rotor degrees 
of freedom for each blade degree of freedom: collective, differential collec- 

tive, lateral cyclic and longitudinal cyclic. In hover, the collective and 
differential collective degrees of freedom are linearly independent of each 
other and of the two cyclic degrees of freedom. The two cyclic degrees of 
freedom are, however, coupled. 

Using quasi-steady or Theodorsen aerodynamics for the fixed hub, multi- 
bladed rotor considered in this investigation, the transformation to multi- 
blade coordinates provides little more information than would an isolated 
blade analysis. The results show the damping of all rotor modes to be the 
same for a single blade mode, and only the frequencies identify the high 
frequency (Q + to) and low frequency (Q - to) cyclic modes. This situation 

i 

occurs because fixing the hub uncouples the blades from each other, and 
because the blades are aerodynamically independent. In addition, note that 
if the modal frequency, o), is greater than the rotor speed, f2; the low 
frequency cyclic mode is the progressing mode and the high frequency cyclic 
mode is the regressing mode. If to < J2, however, the low frequency mode becomes 
the regressing mode and the high frequency mode becomes the progressing mode. 

The nature of the analysis changes when Loewy aerodyiiamics are introduced, 
since the blades are aerodynamically coupled through the wake structure. 

Each of the four independent rotor modes (collective, differential, low 
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frequency cyclic, and high frequency cyclic) resulting from the eigenanalysis 
have different wake geometries associated with them and, therefore, h-va 
different values of the interblade phasing, i|^ . 


ijj = 0 

q 

■ ’’I 

i|j = 2irq/N 

= TTq(l + 2/N) 


Collective 

Differential 

Low Frequency Cyclic 

High Frequency Cyclic 


Since the collective and differential degrees of freedom are independent in 
multi-blade coordinates, there are no difficulties in providing each set of 
equations with the proper aerodynamic coefficients. However, the fact that 
the cyclic degrees of freedom are coupled does present problems. 

If the wake geometries associated with cyclic degrees of freedom are 
used to compute the aerodynamic coefficients for the coupled cyclic equations, 
the same damping is obtained for both the high and low frequency cyclic modes, 
irrespective of the differences in their wake geometries. The reason for this 
occurrence is that in both cyclic degrees of freedom the oscillations of each 
blade lead the one following it by 2TTq/N. Thus, the wake geometries of both 
cyclic degrees of freedom are the same, as are the aerodynamic coefficients. 
Since the hub is fixed, the damping must be the same for both modes. 

There are two ways of circumventing these difficulties. The first is to 
consider only an isolated blade which is coupled to the others through the 
wake geometry [9]. This method is, however, only valid for a rotor with a 
fixed hub. The other way is to devise another transformation for the cyclic 
degrees of freedom which l tcouples them and results in the proper wake geo- 
metries. For this case of the fixed hub, it is immaterial which method is 
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used. The flap- lag perturbation equations can be solved N times Using the 
aerodynamic coefficients for a different mode each time; or be transformed 
into a set of 2N equations, each with the appropriate coefficients, and solved 
simultaneously. 


Results and Discussion 

A sampling of the results obtained from the above analysis are shown in 
Figures 3-7. In these figures, the configuration parameters are as given in 
Table 1, except where noted. 

Figure 3 illustrates the effect of unsteady aerodynamics on a system with 
a low level of damping (note the scale of the ordinate) . While the results 
using Theodorsen aerodynamics do differ from the quasi-steady results through- 
out the frequency range, the differences are small in comparison to the ones 
produced by using Loewy's aerodynamics. The addition of the wake sheets below 
the plane of the rotor appears to have a significant effect on the behavior of 
an isolated blade. There is a sharp spike at the 1/rev lead-lag frequency, an 
N-wave at 2/rev, and smaller dips at 3/rev and 4/rev. 

The results for a four-bladed rote,. (Fig. 4) snow the same characteristics 
as those for the isolated blade. However, the peaks and valleys do not appear 
at every integer frequency for every mode. Instead, the collective mode has 
the 4/rev dip, the differential mode has the 2/rev N-wave, the low frequency 
cyclic mode, has the 1/rev spike, and the high frequency cyclic mode has the 
3/rev dip. Except at the points noted, these curves follow the Theodorsen 
curve quite closely. From other results, it appears that this pattern repeats 
itself at higher frequencies. That is, the collective mode would exhibit 
significant deviations from the Theodorsen curves at 4/rev, 8/rev, 12/rev, and 
so forth. The other modes would behave similarly. In addition, the sequence 
of low frequency cyclic at 1/rev, differential at 2/rev, high frequency cyclic 
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f 


at 3/rev, and collective at 4/rev appears to hold for all other four-blade 
configurations . 

With a system having a low level of damping, like the one above, a 
question arises as to whether the effect of Loewy aerodynamics is significant 
only for that type of system, or for all systems. To check out this possibility, 
a configuration with full elastic coupling and negative pitch-xag coupling was 
analyzed. The results (Fig. 5) show the wake effects to be significant only 
below lead-lag frequencies of 2.5/rev. It appears, therefore, that the woke 
effects become less significant as the level of damping in the system increases. 

In order to determine if a torsionally soft blade would be more sensitive 
to Loewy aerodynamics than a torsionally rigid blade, a pseudo-torsional degree 
of freedom was added to the analysis. While the above analysis does not inclu iu 
the torsional degree of freedom, it can be approximated by using the pitch-flap 
and pitch-lag coupling parameters to produce pitch angle perturbations [11]. 
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^0 


“ 3/'^^0 ^0 


(23) 

(24) 


This approach neglects inertial torsion moments and is generally valid for 
torsion natural frequencies of 5.0 and above. Figure ,6 shows the Loewy curve 
to have much larger deviations from the Theodorsen curve than in any other 
case. It should be noted, however, that the damping values in this case are 
of such a size as to make the assumption of simple harmonic motion somewhat 
questionable. 

Finally, in order to see how the quasi-steady, Theodorsen and Loewy 
results compare for purely simple harmonic motion, a stability boundary was 
plotted for each (Fig. 7). For pitch angles below 0.1 radian and lead-lag 
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frequencies below 1/rev, all three are virtually identical. At higher pitch 
angles, the boundaries of the quasi-steady and unsteady theories diverge from 
one another. The main differences between the Loewy and Theodorsen results 
occur at and above 1/rev. At 1/rev, there is a small dip in the Loewy curve. 
Above that frequency, the curves intersect twice, indicating the type of 
behavior observed in the other figures. 

Another significant feature of Figure 7 is the manner in which the Loewy 
effects do not diminish at the larger pitch angles. It is generally assumed 
that the effects of the wake sheets are more important at small pitch angles, 
where the wake sheets are more closely spaced. In the case of wake-induced 
bending-torsion flutter, decreases in wake spacing are destabilizing [9]. 
However, the pitch-flap problem is linear, and the aerodynamic coupling is 
independent of pitch angle. The flap-lag problem is, on the other hand, 
nonlinear and sensitive to aerodynamic coupling. As the flap-lag aerodynamic 
coupling increases with pitch angle, it overcomes the decrease in the wake 
effect caused by the increased spacing at higher pitch angles, and dominates 
the results. 


Conclusions 

The results of this investigation show that the choice of an aerodynamic 
theory can have a significant effect on stability calculations. Most notable 
is that the addition of wake sheets below the plane of the rotor influences 
the r.otor stability characteristics. While this effect is not exceedingly 
large for .most configurations, except at 1/rev, it does appear that it could 
be quite significant for systems having a low level of damping and for 
torsionally soft blades. 

Secondly, it was shown that the method of multi-blade coordinates is not 
entirely applicable when Loewy aerodynamics are being used. It would also 
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appear that any aerodynamic theory with coefficients that depend on the wake 
geometry would not be compatible with multi-blade coordinates. 

Finally, while the Loewy results are very Interesting, the significant 
effects occur only at integer frequencies. Since rotor systems are designed 
to avoid these frequencies, the wake effects are observed where they have 
minimum influence. Consequently, for most practical configurations, they will 
probably not be quantitatively significant. 
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APPENDIX 

Flap-Lag Equations with Unsteady Aerodynamics 
Equilibrium Equations 
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Table 1. Blade Configuration Parameters 
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Figure 1. — Blade model- 





Figure 3. — Isolated blade (p * ^ lO) at 0^ * 0.1 rad 






Figure 4. — Four-bladed rotor (p = 1.10) at 9^ “ 0.1 rad 










Isolated blade (p - 1.10) at 
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Isolated blade (p = 1.10) with pseudo-torsional flexibility at 
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